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We investigate theoretically the coupling of a cavity mode to a continuous distribution of emitters. 
We discuss the influence of the emitters inhomogeneous broadening on the existence and on the 
coherence properties of the polaritonic peaks. We flnd that their coherence depends crucially on the 
shape of the distribution and not only on its width. Under certain conditions the coupling to the 
cavity protects the polaritonic states from inhomogeneous broadening, resulting in a longer storage 
time for a quantum memory based on emitters ensembles. When two different ensembles of emitters 
are coupled to the resonator, they support a peculiar collective dark state, also very attractive for 
the storage of quantum information. 

PACS numbers: 42.50.Pq, 42.50.Ct, 42.50.Gy, 42.65. Hw 



I. INTRODUCTION 

Understanding the coupling between a cavity and an 
ensemble of emitters was motivated in the early eighties 
by seminal demonstrations of cavity quantum electrody- 
namics (QED) effects jT|. First performed with atoms, 
these experiments were further developed in solid state 
systems, starting with few semiconductor quantum wells 
coupled to planar cavities P]- The interest for this topic 
has been renewed in the framework of quantum infor- 
mation, with proposals to use collections of emitters as 
quantum memories for individual excitations. Indeed, 
ensembles of microscopic degrees of freedom benefit from 
the collective enhancement of the interaction strength [1] , 
while possibly keeping the relaxation properties of a sin- 
gle emitter |3j . This lead to a series of recent proposals 
where cold atoms [1], polar molecules [S] or electronic 
spins |3l|6] coupled to a superconducting cavity have been 
suggested as long-storage quantum memories and optical 
interfaces. This problem also bears some analogy to the 
situation where a nuclear spin ensemble is coupled to a 
single electronic spin [7]. Following these proposals, re- 
cent experiments have demonstrated the strong coupling 
of a resonator to a collection of electronic spins in a crys- 
tal [HE]. However, inhomogeneous broadening is always 
present in the solid state and may eventually limit the 
performance of such a quantum memory. 

In this paper, we study theoretically a cavity coupled 
to a continuous distribution of inhomogeneously broad- 
ened emitters, in the low excitation regime. In the ideal 
case where all the emitters have the same frequency. 
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strong light-matter coupling leads to the formation of 
two polaritonic modes separated by the so-called vacuum 
Rabi splitting [10]. In the situation we aim to describe, 
the emitters bare frequencies are spread over a range that 
can be larger than the cavity linewidth. Our goal is to 
clarify the effect of inhomogeneous broadening on the for- 
mer simple picture in the ideal case, building on an early 
work by Houdre et al In the presence of inhomo- 

geneous broadening, we also find polaritonic peaks. Sur- 
prisingly, their relaxation properties are not only affected 
by the width of the emitters distribution but also by its 
shape. We derive explicit formulas for the polaritonic 
linewidths, showing in particular that provided the spec- 
tral density of emitters in the wings of the distribution 
decays faster than a Lorentzian, the spectral width will 
be dominated by the emitters homogeneous linewidth. 
We call this effect cavity protection. We solve exactly 
the dynamics of the coupled system, showing that, in 
this regime, the two polariton states are well decoupled 
from the other emitters states. As a consequence, cav- 
ity protection reduces very significantly the relaxation of 
an excitation, when stored in one of the polariton states, 
which opens a promising path towards solid-state quan- 
tum memories [T2J US] ■ We finally propose another po- 
tential application of cavity protection, by considering 
a cavity coupled to two inhomogeneously broadened en- 
sembles of emitters. Indeed, this system supports a col- 
lective dark state, which is particularly attractive for the 
storage of quantum information. 

The paper is organized as follows. In Sec|ll]we present 
the model leading to Heisenberg equations in the low- 
excitation regime and obtain an expression for the com- 
plex transmission of the cavity. This expression is an- 
alyzed in detail in Sec |III| where we explore criteria for 
strong-coupling regime, taking into account inhomoge- 
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neous broadening. The transmission pattern allows to 
introduce the notion of cavity protection, whose physi- 
cal origin is analyzed from two different perspectives in 
Sec |IV| and \V\ Finally in Sec |VI[ we study the potential 
of cavity protection in the framework of quantum memo- 
ries. In particular, we discuss the possibility of exploiting 
a collective dark state to store and retrieve quantum in- 
formation with high fidelity. 



II. MODEL 



The system under study is pictured in fig{T] It consists 
in a cavity mode a of frequency ujq , that we shall define as 
the origin of frequencies, linearly coupled with a strength 
gk to a distribution of N two-level systems of frequencies 
LOk and damping rates 7. In the regime where the number 
of excitations is small compared to the total number of 
emitters, each two-level emitter is properly modeled by 
a bosonic mode bk (Holstein-Primakoff approximation). 
The total Hamiltonian writes H = Hcav + Hem + Hint, 
with Hcav = hujQa^a, Hem = J2k f^kb\bk and Hint = 
'iJiJ2k9k{a^bk - bla). 




Port 1 
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Figure 1: Scheme of the emitters-cavity coupled system. The 
cavity frequency is luq. The cavity mode is coupled to the 
outside world via two ports labelled 1 and 2, the fc-th two- 
level system has frequency ujk and interacts with the cavity 
mode with coupling constant gk- 



Using well-known input-output formalism [14', we de- 
fine the external fields Ci„ (injected or pumping field), 
Cr (reflected field) and Ct (transmitted field) that lead to 
the damping k of the intra-cavity field. We also consider 
atomic losses 7, i.e. atomic emission in modes other than 
the cavity mode. The Heisenberg equations are written 
in the frame rotating at the frequency uj of the probe, 
yielding 



= - {k/2 + i{uJo - t^)) a - ^/kJ^Ci^ + ^fffc^fc + /a(i) 



6fc = - (7/2 + ^(cjfc - a;)) Sfc - g^a -f /fe(t) 



(1) 



where /a(t), /fc(i) are noise operators allowing the 
preservation of the commutation relations. From this set 
of equations, and as demonstrated in App|Aj it comes out 
that the evolution of the system can be modelled with 
a generalized Hamiltonian -ffe// involving the respective 
complex emitters and cavity frequencies 10^=10^ — 17/2 
and Coq = loq — in/ 2. Consequently, the system made of 
N atoms coupled to a cavity appears to be equivalent 
to an ensemble of -t- 1 coupled leaky cavities, and the 
problem reduces to the study of the classical evolution of 
the field in each cavity. This exact analogy is the basis of 
the model. Taking the average value and solving analyt- 
ically the set of equations in the steady state regime, we 
get the following expression for the complex transmission 
of the cavity: 



-K/2i 



(ct) ^ 

(c„i) ^o-^^-J2k alli^k - w) 



(2) 



We are interested in the very large number of emitters 
iV, so we describe the emitters as a continuous distribu- 
tion with spectral density p'{uj') spread around its cen- 
tral frequency Wc and normalized to 1. The full width 
at half maximum (FWHM) is denoted A, and is used to 
parametrize each distribution. We replace gk by g{uj',x) 
to account for the fact that the emitters at frequency w' 
can have a spread in the coupling constant, so that 



t{iu) = 



-K/2i 



Uq — iK/2 — UJ — N J duj' - 



p'("') 

■-i'r/2- 



J g^{uj' , x)dx 
(3) 



The integration in x can be carried out independently, 
yielding 



ffo(^') / g^{^',x)dx. 



(4) 



Introducing the collective coupling constant f2 such as 
^l'^ p{uj) = N p' {ui)g'^{uj) with p{u)) normalized to 1, we 
obtain 



K.l2i 



a; - cjQ + in/2 - W{uj) 



with 



W{u) = Vt^ 



p{uj')dijj' 
iu — iu' -\- i^/2 



(5) 



(6) 
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In the following we consider three different continua, 
namely a Gaussian, a Lorentzian, and a rectangular dis- 
tribution. Gaussian broadening is quite common in na- 
ture, from Doppler-broadened lines in gases to e.g. size 
distributions in ensembles of semiconductor nanocrystals 
[l5] and self-assembled quantum dots [16]. Lorentzian 
distributions can be found in certain solid-state systems, 
such as spin ensembles in dipolar interaction |17| or di- 
lute optically active impurities in crystals [TH]. Finally, 
the rectangular distribution is a prototypical example of 
finite bandwidth distribution. The results obtained in 
this case can for instance qualitatively be applied to di- 
lute ensembles of fluorescent molecules in organic crys- 
tals [1^ . For these three distributions, we have obtained 
analytical expressions for the function W{ij), which are 
detailed in Appendix [B} 



poral domain. A first rough interpretation is that strong 
coupling is reached when dephasing processes, that take 
place on a timescale A^^, are slower than than energy 
exchanges, whose period still scales like Note that 

the Rabi period is a collective quantity involving all the 
emitters, even emitters which are not spectrally matched 
to the cavity mode. This apparently puzzling feature had 
already been evidenced in [TT] and is due to the fact that 
the mode interacts with a collective state of the matter 
field. 
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III. PROPERTIES OF THE TRANSMISSION 
FUNCTION 

In this section we discuss the properties of the trans- 
mission function (eq.([5|), in the resonant case. First we 
recall some well-known results in the absence of inhomo- 
geneous broadening (A = 0). In that case, the distribu- 
tion p{uj) is well described by a Dirac delta function, lead- 
ing to W{ll!) = rt^/ {uj + ij/2), and the transmission func- 
tion has two poles X± = zt-y/fi^ _ ((^ — 7)/4)2 -|- «— 

pUj . Strong coupling is reached if >• k,7 and 
is manifested by the appearance of a doublet in the 
transmission pattern located at ±17 (at first order in 
K/n,j/n). These two peaks are the spectral counter- 
part of the coherent and reversible exchange of a quan- 
tum of energy between the cavity field and the sym- 
metrical state \S) of the emitters ensemble, defined as 
\S) = f^^^ X] ffe^l- The transmission coefhcient 
is proportional to the Fourier-Laplace transform of the 
field's amplitude in the cavity initially fed with a single 
excitation {l,G\e~^"'^ff*^'^\l,G) (this result is demon- 
strated in Appendix [A] generalizing ref . j21j and is also 
valid in the case where A > 0). The so-called collective 
Rabi oscillation takes place at the frequency fl defined 
above, which in that case simply equals fl = g^^/N, and 
is damped on a timescale given by the finite linewidth of 
the peaks. In that temporal picture, strong coupling is 
reached when the excitation is exchanged several times 
before being lost in the environment. 

We now study how the strong coupling features are 
modified by inhomogeneous broadening. We have plotted 
the transmission in energy |f(aj)p for fl/A ranging from 
to 3.5 in fig(2j To be only sensitive to the influence of 
inhomogeneous broadening, we have kept k and 7 negligi- 
ble with respect to fl. We have considered the three types 
of distributions introduced in Sec. |ll] namely Lorentzian 
(a), Gaussian (6) and rectangular (c). Whatever the dis- 
tribution, two peaks appear in the transmission pattern 
when r2 > A, a signature of Rabi oscillation in the tem- 
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Figure 2: (Color online). Transmission of a cavity reso- 
nantly coupled to a broad distribution of emitters. (a,d) ; 
Lorentzian; (b,e) : Gaussian; (c,f) : Rectangular. We took 
A = 1 MHz, K = 0.1 MHz, 7 = 10"^ MHz. Bottom : = 3.5 
MHz. These values are typical of NV centers coupled to a 
superconducting resonator. 



Inhomogeneous broadening not only states a novel con- 
dition to fulfill to ensure strong light-matter coupling. As 
it eventually accelerates the damping of Rabi oscillations, 
it also leads to the broadening of the polaritonic peaks, 
as it clearly appears in figj2j In particular, the shape of 
the emitters distribution has a dramatic influence. An 
analytical expression for this width can be derived, in 
perturbation with respect to the small parameter A/Q : 
namely, departing from the strong coupling case in the 
absence of inhomogeneous broadening, we evaluate how 
the poles of the transmission function are modified when 
< A <C ri. For the sake of simplicity we consider 
the limit 7 = 0. The case of finite 7 is studied in Ap- 
pendix [C] in the limit 7 ^ A, which corresponds to the 
experimental situations we aim to describe. Using the 
Sokhatsky-Weierstrass formula in eq.([6]) we have 



W{uj) 



p 



p[L0')dljj' 



(7) 



The modified poles of the transmission function are 
expected in the vicinity of ±f2, so that we develop the 
expression ofW{uj) for w ~ 17 >• A : 



(8) 
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yielding for the poles of the transmission function (at 
first order in k/VI and second order in A/f2), A± — 
±J7 + ^ >v+27rO p{n) ^ pij^a^iiy^ keeping a finite 7 leads to 
the modified expression for the full width at half maxi- 
mum of the peaks : 



r= (K + 7 + 27rp(f7)172)/2. 



(9) 



Looking at eq.(|9|, it appears that in the strong cou- 
pling regime, the polaritonic peaks remain located at ±17, 
but that inhomogeneous broadening adds a contribution 
to their linewidth. This contribution writes 2t:VI^p{Q?) 
and scales like the density of emitters at the real fre- 
quency of the poles. This feature explains the sensitivity 
to the distribution shape that clearly appears in fig(2] 
The polaritonic linewidth decreases upon increasing fJ, 
provided the distribution p{uj) decays faster than 
The Lorentzian distribution is the limiting case for which 
the linewidth tends towards a constant A : whatever the 
coupling, the polaritonic linewidth is governed by inho- 
mogeneous broadening. On the contrary, in the Gaussian 
and rectangular cases, increasing the ratio J7/A allows to 
get rid of the influence of the parameter A, so that the 
width of the peaks only depends on the losses of the cav- 
ity and of individual emitters. In the rectangular case, 
this ideal behavior is even reached for finite values of the 
collective coupling strength VL (while it remains a limit 
in the Gaussian case). This effect, that we call cavity 
protection, leads to an enhanced lifetime of the Rabi os- 
cillation and has interesting consequences for quantum 
information storage as we show in SectiorjVlj 



IV. ORIGIN OF PEAK BROADENING 

Before focusing on applications opened by cavity pro- 
tection, we give an interpretation of peaks broadening. 
This amounts to understanding the damping of Rabi 
oscillations, which occurs even in the absence of any 
radiative losses k = 7 = 0. Our approach is based 
on a seminal paper of Fano [22j, and consists in the 
diagonalization of the total Hamiltonian of the system 

In the absence of inhomogeneous broadening, prepar- 
ing the system in the initial state |1,C?) gives rise to 
Rabi oscillations between the atoms and the field. This 
state is a coherent superposition of two eigenstates of the 
Hamiltonian, namely the polaritons |V'±) — ^ |0, S") ± 

|1,G), of energies ±fi.r2, where \S) is the symmetri- 
cal matter state defined in Section III. Rabi oscillation 
is a quantum beat between these two components. In 
particular, all other emitters states, which do not inter- 
act with the electromagnetic field and are usually called 
"dark states", remain uncoupled. The presence of inho- 
mogeneous broadening strongly modifies the features of 
the emitters cavity coupling. Introducing the continuous 



basis of bare emitters states \uj) of energy huj, we write 
the matrix elements of H as 



{l,G\H\l,G)=h^o 
{J\H\1,G) = hVl^p{Lo') 
{uj'\ H jcj) = hS{u! — U!')UJ , 



(10) 



where the coupling is normalized per unit frequency. An 
eigenvector \ipu) of H with energy huj is searched under 
the form 

= a{uj) |1, G) + / doj'b{u, c^') , (11) 



where the quantity |a(a;)p is normalized with respect 
to Lu. For distributions whose support is not bounded, as 
it is the case for Lorentzian and Gaussian, the solution 
of the eigenvalue equation has been carried out by Fano 
in [22 , yielding the normalized eigenvectors: 



v/pRf7 (^|l, G)+Pf 1^')) + Ciu) lu) 



^G{ojy + (7rp(w)r!2)2 



where Ff stands for principal value and 



(12) 



(13) 



The amplitude of probability to find the excitation in 
the cavity mode can finally be written 

(l,G|e-*^*/''|l,G) = (l,G|e-'^*/'"' / rfc^'a*(w') IV"-') 



da;'|a(tj')|^ 



(14) 



It can easily be shown that |a(w)p is proportional to 
the transmission coefficient in energy |i(w)P (namely, 

|a(w)|2 = 172p(w) for 7,K ^ 0), so that |t(cj)|2 

corresponds to the Fourier transform of the occupation 
amplitude of the cavity mode. As we have checked in 
Appendix |D] this result is completely consistent with the 
formalism of Laplace transform used in Section III in the 
absence of external sources of losses. 

This approach sheds new light on the transmission 
function studied in Section III, which directly reflects 
the overlap between the initial state |1, G) and the con- 
tinuum of eigenstates jV'i^) of the Hamiltonian. The 
two peaks characteristics of the strong coupling regime 
show that this initial state is a coherent superposition of 
two wave packets, reminiscent of the polaritons obtained 
when A = 0. As the eigenstates of the Hamiltonian form 
an infinite continuum, these wavepackets always have a 



finite width, responsible for the damping of Rabi oscilla- 
tions. Nevertheless, as it was shown above, increasing the 
collective coupling may drastically change the shape of 
this overlap and eventually, lead to the narrowing of the 
peaks for distributions p(ui) decaying faster than a 
phenomenon that was defined above as cavity protection. 

Distributions with a bounded support of width A (rect- 
angular for example) , provide an interesting limiting case 
where cavity protection is almost perfect. As a mat- 
ter of fact, if > A, the Hamiltonian eigenstates not 
only consist in a continuum lying within the sup- 
port of the distribution, but also in two discrete states 
IV'-i-) and IV'-), located around uj = if] (at first order 
in A/f2), corresponding to the polaritons ]'(/'+) and 
when A = 0. The initial state 11,0) mostly overlaps 
with these two eigenstates, making the problem similar 
to the case of standard Rabi oscillations in the absence 
of inhomogeneous broadening. In particular, if p{uj) is 
rectangular, the overlap of |1, G) with the discrete states 
equals C — 1 — (l/8)(A/r2)^, giving rise to Rabi oscilla- 
tions of infinite duration characterised by a contrast C. 

To conclude this part, we emphasize that the total 
damping rate F — (k + 7 + 2tt p{Vl)Vl^) / 2 evidenced in 
Section III shows contributions of essentially different na- 
ture. The first type, related to n and 7, is due to the 
irreversible loss of the excitation in the environment of 
the cavity or the emitters. The second type, related to 
7rp(f2)f2^, is Hamiltonian and thus reversible in princi- 
ple with CRIB experiments. It is due to the interaction 
of the cavity with a continuum of emitters, leading to 
progressive dephasing of Rabi oscillations. 



V. OPEN SYSTEM APPROACH 

The approach developed in Section IV gives an inter- 
pretation of the peaks broadening within a Hamiltonian 
formalism. In this part, we adopt another point of view 
based on quantum open systems. As it was exposed 
above and pictured in figjSj in the absence of inhomoge- 
neous broadening, the symmetrical state | S) is decoupled 
from the dark states. The excitation initially injected in 
the cavity mode remains thus trapped in the " small sys- 
tem" consisting in the two polaritons | ')/'+) and IV"")- 
When inhomogeneous broadening is switched on, the 
symmetrical state couples to the dark states, which ap- 
pear as an environment in which the excitation can decay. 
Broadening of the polaritonic peaks can be attributed to 
the decoherence induced by the bath of dark states. This 
picture is inforced by the computed expression for the 
width of the transmission peaks, F = 27rri^p(il), which 
could be interpreted as a natural linewidth for polaritons 
"dressed" by the environment of dark states. Neverthe- 
less, the analogy should be used with caution, as the 
coupling with the bath is not Markovian. This naive pic- 
ture has still the advantage to give an intuitive insight 
on cavity protection, which is nothing but energetically 
decoupling the polaritons from the bath of dark states. 
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as initially suggested in [TB] . 




|0,G) %G) 

Figure 3: (Color online). Schematic diagram of the open 
system approach to inhomogeneous broadening, (a): A = 0, 
the states IV'i) a-re isolated from the degenerate dark states 
\uS). (b): A ^ 0, the states IVi) are coupled to the \uj) 
states, which are non-degenerate in this case, with a coupling 
strength proportional to A . 
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Time (units ot l/A) Time (units of 1/A) Time (uuils of 1/A) 

Figure 4: (Color online). Probability to recover an excitation 
initially stored in the state \'4>+) after a time t. We took A = 1 
MHz, K = 7 = 0. (a) : Lorentzian ; (b) : Gaussian ; (c) : 
Rectangular. Red Dashed line : Q, = 1 MHz; Green Dotted 
line : Vl^2 MHz; Blue line : Q, = A MHz. 

To study the dynamics of the polaritonic relaxation, 
we have exactly computed the evolution of the state of 
the system initially prepared in for different val- 

ues of the collective coupling strength Vl, and for the 
three types of distribution, keeping the same FWHM 
A = IMHz. We have plotted in fig|4] the probability 
I ("0+1 e~'^^'^i^ P of finding the excitation in the po- 
lariton, as a function of time. For the sake of clarity, we 
have neglected again the losses k = 7 = (realistic val- 
ues are considered below). As it can be seen in the figure, 
if the distribution is Lorentzian, the excitation exponen- 
tially decays in the environment, whatever the coupling 
51, on a typical timescale A~^. This is consistent with the 
spectral study performed in Section HI, where the width 
of the polaritonic peaks does not depend on the coupling 
with the cavity. On the contrary, the effect of cavity 
protection can be observed on the two other distribu- 
tions. Damping is strongly inhibited as soon as f2 > A if 
the distribution is Gaussian, but is always present what- 
ever the coupling, which is the counterpart of the finite 
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linewidth of the transmission peaks. Finally, in the case 
of a rectangular distribution, two timescales are visible. 
The initial state mostly overlaps with the discrete 

state \ip+) defined above, but also with the continuum 
of eigenstates {tpuj). The coherent superposition of the 
continuum of frequencies is damped on a short timescale 
A~^, so that the probability quickly converges towards 
the quantity | ^V+l which also scales like (A/O)^. 

VI. APPLICATION TO QUANTUM MEMORIES 

The previous Sections establish that for distributions 
allowing cavity protection, increasing the collective cou- 
pling f2 dramatically increases the potential storage time 
of one excitation in the polaritonic states, as energetic 
decoupling from the dark states is more pronounced. In 
particular, this storage time becomes insensitive to de- 
phasing processes induced by inhomogeneous broaden- 
ing. This allows to treat an inhomogeneous distribution 
as an effective oscillator of ground state \G) and first ex- 
cited state 15*), that benefits from the collective coupling 
il to the cavity and whose relaxation properties are solely 
governed by individual emitter properties 7. As a con- 
sequence, cavity protection opens the path to the imple- 
mentation of long lived solid-state quantum memories, by 
exploiting ensembles of microscopic degrees of freedom, 
whose coherence times are remarkable. In this Section 
we use our modelling to estimate the performances of 
two such types of quantum memories. 

A. Quantum memory based on dispersive coupling 

Here we evaluate the potential of a broad ensemble of 
emitters dressed by a cavity mode for quantum informa- 
tion storage. The coupling should be dispersive to freeze 
Rabi oscillations between the mode and the atoms. This 
system offers an interesting situation where information 
has to be protected against two types of losses : the cav- 
ity losses, which are more critical when the mode and 
the distribution of emitters are on resonance, and the 
losses in the dark states, which on the contrary, are all 
the weaker as the atoms-cavity detuning is smaller. The 
atoms-cavity detuning is thus the result of a tradeoff, and 
can be optimized with our modelling, as we show below. 

The protocol of the quantum memory is the following. 
First, the detuning d between the mode and the center 
of the distribution is slowly swept from —00 to a finite 
positive value, thus adiabatically mapping the quantum 
state of the cavity mode onto the emitter's ensemble : 
{a |0) + /3 |1)) |G) ^ |0) (a \G) + p jV'+l^))). We have 
introduced the dressed state = cos(0/2) |0, 5) + 

isin(0/2) |1,G), and the mixing angle cot(0) — (5/(2Sl). 
The transfer of the excitation should be realized on a 
timescale longer than the Rabi period, but shorter than 
so that no dephasing mechanism affects the pro- 
cess, this can be achieved under strong coupling as in 



this case ^ A . The expected fidelity F{t) of such 
a quantum memory can be exactly computed with the 
present model; in particular, in the case where a single 
photon state is stored (/3 = 1), we get the simple expres- 
sion J" = I {ip\{8)\ e-^"-ff* p. We have plotted 
this quantity in figjsj As explained above, must be 
optimized by properly choosing the detuning S, which 
should be low enough to maintain cavity protection, and 
high enough to reduce the sensitivity to cavity losses, 
which typically scale like K;(r2/(5)^. The maximal detun- 
ing Sm leading to an efficient protective energy gap is 
H'^/Sm ^ A [12]. This condition induces an optimal re- 
duction of the cavity losses by a factor of (fi/A)^. 

The trade-off in the detuning clearly appears in the in- 
set of figjSj where we have plotted J^, as a function of the 
detuning S, after ten cavity lifetimes, for different values 
of the ratio fl/A. We have used standard parameters for 
circuit QED technology [8]. As it appears in the figure, 
a quantum memory based on a Gaussian distribution of 
emitters of linewidth A = 1 MHz, strongly coupled to 
a cavity of width «; = 0.1 MHz with a strength ft = 40 
MHz would yield a typical fidelity of 90% after 100 ^s, a 
remarkable storage time compared to the lifetime of the 
cavity mode (10 fis) and the typical dephasing time of 
the ensemble (1 fis). 
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Figure 5: (Color online). Maximized fidelity of regaining 
the excitation initially stored in the state \tp+) after r = 10 
cavity lifetimes, as a function of Q/A. We took A = 1 MHz, 
K = 0.1 MHz, 7 = 10~* MHz. Inset : same quantity £ as a 
function of detuning 5, after t. Green dotted line : Q — 40 
MHz ; Black solid line : fl = 20 MHz; Blue dashed line : 
Q = 10 MHz; Red dotted fine : Q = 5 MHz. 



B. Quantum memory based on two emitters 
distributions 

We focus now on a second type of quantum memory, 
based on two distributions of emitters allowing cavity 
protection, respectively detuned by +6 and —S with re- 
spect to a cavity. Note that the case of a mode cou- 
pled to two such discrete emitters of ground and excited 
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states \gi) and jci) is exactly solvable, the poles of the 
transmission revealing the complex eigenfrequencies of 
the system [53]. In particular, when the emitters are on 
resonance with the mode {5 = 0), the antisymmetrical 
state {\ei,g2) — j^i, e2))/-\/2 is not coupled to the elec- 
tromagnetic field. This dark state is naturally protected 
against spontaneous emission in the cavity, a property 
that can be used to store quantum information during 
a typical timescale given by the atomic dephasing time. 
Note that for artificial atoms like superconducting qubits 
or quantum dots this time can be quite short, which is a 
severe drawback for quantum computation on chip. Here 
we suggest an experiment to prepare and exploit this dark 
state as a quantum memory, in the case where the dis- 
crete emitters are replaced by broad assemblies of atoms. 
This proposal allows to benefit from the collective atoms- 
cavity coupling, while the storage time now corresponds 
to the dephasing time of individual emitters, and is thus 
potentially quite long. Note that this idea is typical of 
the so-called hybrid circuits technology [3]-(B]. 

First we have checked the validity of the effective model 
if two ensembles are coupled to the cavity. We have plot- 
ted in figj6^ the exact transmission |t(cLi)p of a cavity 
coupled to two Gaussian ensembles and verified that the 
position of the peaks are fitted by the eigenenergies com- 
puted in the discrete case. Moreover, we have superim- 
posed the transmission resulting from the exact calcula- 
tion and from the discrete model, as it can be seen in 
figj6]3 after focusing on the central peak of the transmis- 
sion pattern: the excellent agreement between the two 
plots fully validates the effective approach. This central 
peak corresponds to the eigenstate \tl>d) resulting from the 
coupling between the cavity mode and the antisymmet- 
ric state 1^) = (|Gi,S'2) — jS*!, G2))/\/2, its expression 
being l^d) = (z,5|l,Gi,G2) + r!^|0,^))/VFTW. 
When (5 3> ri, the excitation is mostly in the cavity, and 
mostly in the matter field in the opposite case. This 
change of nature clearly appears in the narrowing of 
the peak while lowering (5, as it can be seen in the fig- 
ure, and confirmed by the expression for its linewidth 
Td = + 2f]2^)/(52 + 2f]2). Note that this mod- 
elling might explain some recent experimental results 
[5], in which a superconducting cavity is strongly cou- 
pled to a inhomogeneous ensemble of NV centers of spin 
1. Because of the geometrical strain, the transitions 
|ws = 0) — )■ \ms = 1) and |ms = 0) — >■ \ms = —1) are 
splitted, which can be modeled by two ensembles of emit- 
ters of different central frequencies. The visible presence 
of a narrow peak at the cavity frequency explains quali- 
tatively the effect discussed above. 

Coming back to the general case of two distinct ensem- 
bles, the state \i}}d) could provide a new type of quantum 
memory as mentioned in the beginning of this Section. 
The protocol consists in feeding the cavity mode with a 
single photon while the ensembles are largely detuned, 
thus preparing the state |l,Gi,G2), then adiabatically 
transferring the excitation to \'ipd) by slowly lowering 5. 
Yet the ensembles cannot be brought to resonance with 



the mode as it would be the case for two discrete emitters. 
As it appears in fig|6j:, the effective model breaks down 
when (5 ^ A. At this point indeed, the distributions 
of emitters start to spectrally overlap with the central 
peak, leading to its broadening. This yields a minimal 
linewidth ^ 7 -I- (A^/2r2^)K, allowing to typically re- 
duce the cavity losses by (f2/A)^. Here again the ratio 
(ri/A)^ appears as a major figure of merit for devices 
based on inhomogeneous ensembles strongly coupled to 
cavities. 





-0.4 -0.2 0.2 0.4 -0.4 -0.2 0.2 0.4 

Frequency (units of A) Frequency (units of A) 



Figure 6: (Color online), (a) : Transmission of a cavity cou- 
pled to two Gaussian distributions of emitters, each detuned 
by +5 and —5 from the cavity frequency, 5 is swept from 
to 8 MHz. We took ^ = IMHz, A = O.lMHz, k = 0.5MHz, 
7 — 10~*MHz. (b) : Focus on the central peak with 5 = 0.5 
MHz. Solid red line : Gaussian profile. Blue dashed line : two 
emitters of homogeneous linewidth 7. (c) : 5 = 0.15 MHz. 



VII. CONCLUSION 

We have shown that if an inhomogeneous distribution 
of emitters is strongly coupled to a cavity, the ensemble 
can be treated as a single effective emitter collectively 
coupled to the mode, whose relaxation is governed by 
single emitter's properties, provided that their spectral 
distribution decreases faster than 1 /lj^ . This effect called 
"cavity protection" offers promising perspectives in the 
framework of quantum information with solid state inte- 
grable devices, in particular regarding the implementa- 
tion of long lived high fidelity quantum memories. These 
results are quite general, and can fruitfully be applied 
to numerous important physical systems, ranging from 
semiconductor emitters coupled to optical cavities, to en- 
sembles of spins in circuit QED. 
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^i('S) = (I7 G*! (s + i-ffe/z/fi) ^|1,G). Inverse Laplace 
transform of this coefficient gives back the quantity 
(1, G\ e~^"'fi/^'^ |1, G). We easily get 



Appendix A: Dynamics 

In this part, we establish the link between the complex 
transmission of the cavity, and the evolution of the sys- 
tem if the mode a is initially fed with a single photon. 
This evolution is governed by the set of equations ([T]) 
written in the free frame (w = 0). The input fields are in 
the vacuum, the state of the system is |1, G) = {0) |0) 
where |0) is the ground state of the total system. We are 
interested in the quantities (a(t)a'^(O)) and {bk{t)a^ {Q)) , 
which represent the probability amplitude of the excita- 
tion in the cavity mode and in each emitter respectively, 
as it will appear later. The average values are taken in 
state |0). We get 

(d(t)at(O)) = - + lUJo) {a{t)a^O)) + ^ 9k{bk{t)a\0)) 

k 

{bkit)a\0)) = - (7/2 + iLOk) {bkit)a\0)) - gfe(a(i)a^(0)) 

(Al) 

Defining the vector \tp) of coordinates 
((a(t)at(0)),...,(6fc(i)at(0)),...), its evolution 

simply follows the Schrodinger like equation 



Heff/h = 



I 



\ 



^0 «5i «.92 



(A2) 



We have used the complex frequencies for the cavity 
Wo and for the emitters w^- defined above. Note that 
these results are in full agreement with the ones ob- 
tained in the Green function formalism by Kurucz et 
al [12]. It appears that the dynamics of the problem 
can be modeled with the effective Hamiltonian i?e/ / . In 
particular, one can define an effective evolution opera- 
tor 0{t) = e»ffe//t/ROe-*^="*/'', such that (a(i)a^(O)) = 
(0| a(0)e-*^=//*/'^a'f(0) |0). This quantity can be rewrit- 
ten (1, G| e-^^-ss^l^ |1, G), justifying that we talk of the 
probability amplitude of the excitation in the cavity 
mode, starting from the initial state |1, G). The problem 
is solved using e.g. standard Laplace transform method. 
Defining £ (/(t)) = F{s) = exp{st)f{t)dt, we have 



his) 



9k 



(A4) 



From eq.([5]) and eq.(A4), we finally write the link be- 
tween the transmission coefficient in amplitude t(uj) and 
the coefficient ti{s) characterizing the dynamics of the 
system. 



t{uj) ^ --ti{-iuj) 



(A5) 



This establishes the relation between the amplitude 
Qi(<) = (1, G| e~^"'"^/^ |1, G) and the transmission tlu) 



2 . X 

K 



(A6) 



One can use the method exposed above to com- 
pute the expression of the probability amplitude 
for a state initially prepared in |V''].((5) ), namely 



{ipl{S)\e-'"''ff*/^ \^PHS)) studied in Sectic 
eral, we can decompose it as 



VI 



In gen- 



{i^i{s)\u,ff{t)\rAS))^ 

= cos2(0/2) (1, G| Ueff |1, G) + sm^{9/2) (0, S\ U^ff |0, S) H 
+ isin(0/2) cos(f?/2) ((0, S\ U.ff |1, G) - (1, G| U,ff |0, S)) 
= cos^{e/2)ai(t) + sm^{e/2)a2{t) + 
+ isui{e/2)cos{e/2) {a^it) - ai{t)) , 

(A7) 



where Ueff{t) 



We need only to calculate the four matrix elements 
ai{t). Defining ti{s) = C{ai{t)) we obtain in the case of 
a continuous distribution. 



, , W(is) , 

t2{s) = ^h{s){s + lUJo) 



his) = his) 



iWiis) 



(A8) 



tiis) = -his) 



Appendix B: W(cj) for specific distributions 



\m) = {{s + ^H,ff/h)-' \m)) , (A3) 

where we have used the Laplace transform property: 
= s\^is)) - IV'(O)). We finally define 



We now evaluate the function for all the specific 

continua analyzed in the paper. This function allows 
the evaluation of the complex transmission using = 
(K/2i)(a; — ajo+i'*/2 — but also appears in other 
formulae. 
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Gaussian 



The Gaussian distribution writes 

^g-(^^l„2)/A=^^ is thus, 



1 vin2 

~i A 



Remembering that 



duj'- 



, A/Vln2 



erfc(— iz) 



(Bl) 



(B2) 



where erfc is the complex complementary error function, 
it comes 



\/ln 2 _/oj + iT/2^x2 

WG{'^) = -i v^e ^A/v^^ erfc 



A/Vln2 

(B3) 



Rectangular 

In the case of a rectangular distribution, the density of 
emitter is p{uj) = \{<d{i^ - A/2) - <d{uj + A/2)) we have 



2fi2 / A 

WR{io) - — ArcTan — 

I A \7 ^ 2«aj 



(B4) 



Lorentzian 



The density of emitter is p(a;) — 



^ thus 



(A/2)2+t^2 



uj + ^7/2 + iA/2 



(B5) 



From the equation above we see that for Lorentz dis- 
tribution we do not achieve cavity protection, i.e. the in- 
homogeneous broadening always contributes as if it were 
homogeneous. 



f 

w{uj) = n^p 



°° p(uj')duj' 



UJ — uj' 

7 ,r Pi^')d^' 



{uj-uj'y 



(C2) 



We are interested in the development of W{uj) near 
the poles of the transmission function in the absence of 
inhomogeneous broadening, namely uj ^ Q. Denoting 
r = uj'/uj, and using the identity '^r'^ — 1/(1 ~ '')j '^6 
find ; 



Wi^)^^-\l + J2>^^^np{uJ) 



oo 



uj^ 2 \ ^-^ UJ 

\ k=l 



(C3) 



fc 



where pk is the fc-th moment of the distribution p{uj) 
about its origin 



fJ-k 



dujp{uj)uj^ . 



(C4) 



Note that this development is only valid if w ^ w', 
which is the case in the present study as w ^ ^ A > 
uj' . From the normalization and considering only sym- 
metric distributions, we have p^ = \ and pi = 0. p2 
gives the first non-zero correction and it is typically pro- 
portional to the square of the FWHM (as an example, 
P2 — A^/(21n2) in the case of a Gaussian distribution). 
To first non-zero correction we have: 



W{uj) = —{! + P2/0J^) TTQ^piuj] 
UJ \ 2 uj^ 



n^ii + P2/uj^) 

UJ + i^/2 



(C5) 



iiT^l^ p{uj) 



where we have used f2 ^ 7. One easily infers the modi- 
fications to the transmission poles induced by inhomoge- 
neous broadening. They are located at 



Appendix C: Development with finite 7 

We start rewrite W(uj) from eq.([6]) as 

, uj'"^ p{uj' + Uj) 



w{uj) = duj 



iiT^^ / duj' 



uj'^ + 7^ uj' 



7 



7r(w'2 -1-72) 



p{uj + uj') . 



(CI) 



The integrands contain products of a function of width 
7 and another with width A. If 7 ^ A, the integrals take 
the form: 



Their width check F 



, m correspon- 



dence with what stated in Sec. |III| Note that this proce- 
dure is only valid for distribution with well defined mo- 
ments. This is not the case of the Lorentzian, neverthe- 
less W{uj) can be exactly evaluated in this case. The 
exact calculations for the 3 cases taken under considera- 
tion are the subject of appendix [BI 
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Appendix D: Two ways to obtain the temporal 
evolution 



We have found two ways to evaluate ai(t) 
(l,G|e-*^="*/''|l,G), the first at Appendix |A| 
uses a Laplace-Fourier transform of ~t{uj)/{K/2 
the second uses the standard Fourier transform o 

as in Sec. 



IV 



2TTn^ p{uj)\t{uj)/{K/2)\^ for K,7 

The first way is more general in the sense that it can 
include emitter and cavity radiative losses, the second 
describe a reversible process originated in a Hamiltonian 
evolution. We now show that both ways coincide when 
we disregard losses. 

From Appendix [A| we have 



ai{t)e''^*dt = ti{-iuj) , 



(Dl) 



where, if 7, k — >■ 



if we consider time reversibility of the lossless dynamics 
we have ai{—t) = {ai{t))* and thus 



(D4) 



ai(t)e*"*dt. 



Eq.(D3l and eq.(D4) together give 



ti{—iuj) 



(D2) 



We now take the real part of eqs.(Dl , D2| , yielding 



Lj-LJo- n^pj ^^^^ + i^fiv(w) 



/oo 
ai{t)e"^'dt = 27rn^ p{uj) \ti\^ , (D5) 
-00 



which is precisely what we find applying the inverse 



Fourier transform in eq.(14). Note we had to use the 



(D3) time-reversibility which is only valid in the lossless case 
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